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ABSTRACT 



_ ■ We present limits to anisotropies in the cosmic microwave background 

<D • radiation (CMB) at angular scales of a few arcminutes. The observations 

■ were made at a frequency of 142 GHz using a 6-element bolometer array 

pq ■ (the Sunyaev-Zel'dovich Infrared Experiment) at the Caltech Submillimeter 

Observatory. Two patches of sky, each approximately 36' x 4' and free of known 
>- ■ sources, were observed for a total of 6-8 hours each, resulting in approximately 

80 independent If 7 full- width half-maximum pixels. Each pixel is observed 
with both a dual-beam and a triple-beam chop, with a sensitivity per pixel 
Q ■ of 90-150 /iK in each chop. These data have been analyzed using maximum 

t^^ ■ likelihood techniques by assuming a gaussian autocorrelation function for 

■^ . the distribution of CMB fluctuations on the sky. We set an upper limit of 

Q_|. AT/T < 2.1 X 10^^ (95% confldence) for a coherence angle to the fluctuations 

Q ■ oi I'.l. These limits are comparable to the best limits obtained from centimeter- 

^ . wavelength observations on similar angular scales but have the advantage that 

cd . the contribution from known point sources is negligible at these frequencies. 

>• . They are the most sensitive millimeter-wavelength limits for coherence angles 

^ . < 3'. The results are also considered in the context of secondary sources of 

j^ . anisotropy, speciflcally the Sunyaev-Zel'dovich effect from galaxy clusters. 
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Subject headings: cosmic microwave background - cosmology: observations 



1. Introduction 

Measurements of anisotropies in the spatial distribution of the cosmic microwave 
background radiation (CMB) are a powerful probe of the early universe. In the standard 
inflationary model, anisotropies were imprinted on the CMB when the universe combined 
at 2; ~ 1000 (for a review of cosmological theories see, for example. White, Scott and 
Silk 1994). The attenuation of these primordial fluctuations in the CMB by the finite 
thickness of the surface of last scattering, and photon diffusion, suppress the anisotropy 
power spectrum at arcminute angular scales. The exact shape that is predicted for the 
power spectrum in this regime is strongly dependent on the assumed cosmological model. 
Additionally, an early period of re-ionization may have only a small effect on the anisotropy 
power spectrum at degree scales while strongly affecting the magnitude and shape of 
the power spectrum on arcminute scales. Thus measurements of, or upper limits on, the 
magnitude of the CMB power spectrum at arcminute scales are a powerful discriminant 
between competing models. Secondary sources of arcminute-scale anisotropies include the 
Sunyaev-Zel'dovich (S-Z) effect — the inverse Compton scattering of CMB photons by hot 
gas residing in the potential wells of galaxy clusters. Measurements at arcminute angular 
scales can provide important information for theories of large-scale structure formation via 
number counts of galaxy clusters with a measurable S-Z effect. 

The Sunyaev-Zel'dovich Infrared Experiment (SuZIE) is a 6-element bolometer array 
that has been used to make the first detections of the S-Z effect at millimeter wavelengths 
( IWilbanks et al. 1994| , Polzapfel et al. 1997a|) . Coupled to the Caltech Submillimeter 



Observatory (CSO), this array has a sensitivity of 0.3 Jy/vHz at 142 GHz in each of six 
pixel pairs and angular resolution of ~ 1'7. This, and a scan strategy designed to minimize 
systematic effects, make it an ideal instrument with which to search for arcminute scale 
anisotropies in the CMB. 

This paper describes observations made with the SuZIE receiver of two regions of sky 
that are free of known sources. The data are used to set limits on CMB fluctuations that are 
assumed to be distributed on the sky with a gaussian autocorrelation function (GACF). A 
companion paper ( |Ganga et al. 1997aD considers the application of these results to specific 
cosmological models. In this paper we also consider the limits that these observations place 
on number counts of S-Z clusters. 
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2. The Instrument 



The SuZIE bolometer array (described in detail in [Holzaptel et al. 1997b| ) is operated 



at the Caltech Submillimeter Observatory (CSO) on Mauna Kea. At 142 GHz, emission 
from the sky and the warm telescope optics are measured by SuZIE to contribute 36 K in 
brightness temperature. A tertiary mirror re-images the Cassegrain focus of the telescope 
to a focal plane comprising two rows of three bolometers cooled to 300 mK. Each bolometer 
is coupled to the primary mirror by a Winston concentrator that defines a If 7 full-width 
half-maximum (FWHM) beam on the sky. The two rows are separated by 2', and adjacent 
pixels within a row by 2'3. The illumination pattern of each Winston cone on the primary 
mirror is controlled by a 2 K Lyot stop located at an image of the primary produced by the 
tertiary mirror. To reduce systematic effects from spillover, only 8 m of the 10.4 m primary 
diameter are used. 

SuZIE can be configured for operation at 142, 217 or 268 GHz by placing an appropriate 
metal-mesh band-pass filter in front of the Winston cones. The observations described 
here were made at 142 GHz with an 11% band-width filter. The spectral response of the 
instrument, in its telescope-ready configuration, was measured in the laboratory prior to the 
observations, including measurements that demonstrate negligible out-of-band leaks. The 
measured scatter of the band centers of all six channels about the mean value is less than 
1.5% and that of the band- widths is less than 5%. The optical efficiency of each channel, 
including the effects of all filters, is measured to be 37%. 

Electronic differencing between pixels in the same row is carried out to remove common- 
mode response to atmospheric and telescope emission. Two differences corresponding to 2'3 
on the sky and one to 4'6 are obtained from each row of the array. This is accomplished by 
placing the two bolometers in an AC-biased bridge circuit ( |Wilbanks et al. 1990| ), the output 
of which is synchronously demodulated to produce a stable DC signal corresponding to the 
brightness difference on the sky. In terms of atmospheric subtraction, this is equivalent to a 
square- wave chop on the sky at infinite frequency. High rejection (> 40 x) of common-mode 
signals from detector temperature fluctuations and telescope and atmospheric emission 
( plezer, Lange fc Wilbanks 1992]) is achieved by trimming the amplitude of the bias 



voltage, and thus the responsivity of each detector. To ensure that the pixels that are being 
differenced are viewing the same column of atmosphere, and thus that maximum rejection 
of common-mode atmospheric fluctuations is obtained, the tertiary mirror is designed to 
maximize the overlap of the beams from the six Winston cones on the primary mirror 
without significantly degrading the focus. The illumination of the primary was measured 
at 217 GHz using a mobile 80 cm^ load, and was verified to fall to < — 60db at a radius of 
4.0 m from the primary center. 
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3. The Observations 

SuZIE is designed to measure the S-Z effect in galaxy ciusters. As part of this program, 
regions of sky free of known sources are observed to provide a check for baseline effects that 
might contaminate the S-Z measurements. We selected two such regions for observation in 
April 1994 at the positions listed in Table ^ using IRAS catalogs and the NRAO 5 GHz 
survey ( Pecker, White and Edwards 199"l|) to avoid known sources. Each field is near a 
SuZIE target cluster, but at least 10 core radii away from the cluster center. The S-Z 
contribution to each pixel from any residual hot gas at this distance is ATcmb < 3/iK, 
assuming an X-ray core radius of 1', a ^/-parameter of 3 x 10~^ and, pessimistically, that the 
cluster gas is described by an isothermal, (3 = 2/3, model out to many core radii (see for 
example [Jones and I'brman 1984| ). Including the effects of differencing pixels in the focal 
plane will reduce this still further, depending on the angular scale on which the cluster gas 
is clumped. 

In order to minimize systematic errors that might otherwise arise from position- 
dependent variations in telescope spillover, SuZIE observations are made by fixing the 
telescope position relative to the earth and allowing the source to drift across the array as 
the sky rotates. Each scan is two minutes long, yielding a 30' x cos 6 strip of sky, where S 
is the source declination. After each scan is completed, the dewar is rotated to keep the 
long axis of the array aligned in the direction of right ascension. The starting positions of 
successive scans alternate between 12' and 18' in RA ahead of the location listed in Table |l|, 
providing a check for systematic effects that depend on position within a scan rather than 
on the pointing location of the telescope. Because there are two rows of detectors in the 
focal plane, the complete data set for a single field covers two strips of sky each 36' cos (5 
in RA by If 7 in declination, with twice as much observing time allotted to the central 
two-thirds of each strip as for the outer portions. The observing scheme is summarized in 
Figure |1|. Because both sources are close to the celestial equator, cos 5 ~ 1 and the total 
area of sky covered is 0.057 square degrees. 



4. Data Reduction and Calibration 

Seven data sets were obtained, each corresponding to one night of observation on 
one field (typically 50-100 scans). Each data set is first split to separate the scans that 
correspond to the different RA offsets. The 6 differential measurements per data set that 
are obtained by electronically differencing each detector with other detectors in the same 
row are denoted by pairs dpq where pq = 12, 23 and 31 for the first row and 45, 56 and 
64 for the second row; the differential signal is the result of subtracting the gth from the 
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pth detector voltage. The data were calibrated using drift-scan observations of Uranus for 
which the expected flux was calculated by convolving the spectral model for the source 
( [Grittin and Urton 199tJ| ) with the instrumental spectral response. The uncertainty in the 



absolute flux of Uranus at 142 GHz obtained by this method is /u = ±6%. Corrections 
for atmospheric opacity due to the difference between the elevation of Uranus and the 
elevation of the Fields 1 and 2 observations were < 2% (based on data from the CSO 225 
GHz r-monitor extrapolated to 142 GHz). The result of the calibration is then a signal 
whose units are flux difference between the two beams. This is converted to ATcmb (the 
thermodynamic temperature difference between the two beams) as follows: 



where x = hu/kTcuB and: 
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The temperature of the CMB is taken to be 2.726 K ([Mather et al. 1994|) ; the effect of 



the 0.01 K uncertainty in Tcmb is negligible compared to other sources of calibration 
uncertainty. 

The solid angle of one pixel, fibeam, is calculated from measurements of the beam proflle 
made using Jupiter and Uranus. The 2-dimensional shape of the beam was determined from 
drift scans across Jupiter, with the array offset in declination from the source (transverse to 
the scan direction) by steps of 15". Jupiter is used for this measurement because it is very 
bright at 142 GHz and can thus be used to obtain high signal/noise measurements of the 
sidelobes and the aspect ratio of the beams. However, its large size (40") compared to the 
SuZIE beam makes it unsuitable for a direct measurement of the beam solid angle. This is 
determined instead from drift scans across Uranus, by flrst assuming a circularly symmetric 
beam on the sky and integrating the beam proflle and then correcting for the measured 
aspect ratio of the beams. Note that the beams are very nearly circular (see the Jupiter 
beam map in pJolzapfel et al. 1997b| ) with an aspect ratio (deflned as FWHM in the scan 



direction divided by FWHM in the cross-scan direction) that varies from 0.94-1.13. 

The uncertainty in the beam shapes, /b, (estimated from the scatter of the 6 measured 
solid angles about their mean value) contributes a further 5% to the calibration uncertainty. 
The total calibration uncertainty is then /c = (/^ + fiY^'^ = 8% (Icr). This uncertainty is 
included in the likelihood analysis of the data carried out in §|^. 

Because each detector contributes to the signal measured by two pairs, the six data 
streams are not independent. To remove the degeneracy, the two 2f3 pairs from each row 
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are differenced to generate a double difference: 

tl23 = {di2-d23)/2 (3) 

t456 = (C?45 - C?56)/2 (4) 

The data then comprise 4 independent differences, two single differences, d^i and d^^, and 
two double differences, ti23 and ^456- 

Cosmic ray impacts on the bolometers cause spikes in the data stream that must be 
removed. A point-by-point differentiation is carried out and spikes are identified by a large 
positive excursion followed by a large negative excursion, or vice versa. Approximately 5% 
of the data are removed by this process. The data are then binned by averaging together 15 
samples (3 s of data), equivalent to a 0f75cos5 portion of the drift scan, thus oversampling 
the lf7 beam full- width half-maximum (FWHM) by a factor of 2.3/cos5. 

A best fit offset and linear drift are removed from each scan prior to co-adding all scans. 
The linear drift arises from common-mode variations in signal, primarily from changes in 
the atmospheric or telescope temperature, that are not completely removed by differencing. 
The average values, and variance, of the removed linear drift for each chop are shown in 
Table |^. The effect of this process on the sensitivity of the system to CMB fluctuations is 
small, but is taken into account in the analysis in §|^ 

Because the residual noise is dominated by atmospheric emission fluctuations with 
a l//-type spectrum ( |Holzapfel et al. 1997b| ), the statistical error obtained at each point 
from the binning process is not a good estimate of the sensitivity of the system over longer 
integration times. The correlation time of the atmospheric noise can be estimated from the 
correlation functions of individual scans and is of order 5 s (defined as the time to the first 
zero in the correlation function). Consequently, the binned data points within a single scan 
are correlated over many bins. The effects of this correlation on the statistical properties 



of the data are considered further in § |6.3| . In order to assign a sensible weight to each 
scan, the rms scatter, aj, of the binned data in the jth scan about the best fit offset and 
linear gradient is taken to be a representative measure of the noise in that scan. If the 
atmospheric noise is normally-distributed within a scan and uncorrelated between scans, 
this is a good assumption. Note that this process also assumes that the contribution to the 
rms from true temperature anisotropies within a single scan is negligible. The implications 



of these various assumptions are discussed in §6.3 



The scans corresponding to a single RA offset observed on a single night are co-added 
by combining all of the data at the ith point in each scan to yield a signal yi with 
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uncertainty ai where: 



Ey^,,/^! 
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(5) 



^' = ^TT7T2 • (6) 



Here aj is the rms of the jth scan after removal of the offset and hnear drift, as described 
above. Note that ctj = a, a constant across the entire co-added scan. Figure ^ shows the 
distribution of the aj values for Field 1 taken with the ^64 difference. Note that all of the 
data are included in the co-add and there is no cutting of data based on high values of scan 
rms, aj. Figure ||b shows the distribution of the quantity Vij/a'j evaluated for a single bin, 
again using all of the Field 1 scans and data corresponding to the d^i difference. It can be 
seen that this quantity is normally-distributed with no extreme values or large wings. This 
indicates that the quantity aj is a good estimate of the variance across the entire data set. 

Co-added data from the two RA offsets are then combined by a weighted average of 
data points that correspond to the same position on the sky. Finally, a weighted average 
of the co-added data sets from each night is carried out. Figure ^ shows the final co-added 
data set for Field 1. The rms uncertainty for each 0f75 bin is 130-220 /zK, corresponding 
to a sensitivity of 3-5 mK s^/^ for each difference. Both of the single and double difference 
data sets are shown, with the appropriate beam response to a point source (measured by 
drift scans across Uranus) indicated in each panel. Approximately 6 hrs of integration time 
(190 scans) were obtained on Field 1 and 9 hrs (277 scans) on Field 2. Note that double 
differencing reduces the error bars associated with the ti23 data points relative to those 
of dsi in Figure ^, as residual atmospheric gradients are being removed from the data. In 
the case of t^^Q and ^64, no such improvement is seen, probably due to poor common mode 
rejection in one of the two pairs that are combined to form t^^Q. 



5. Statistical Analysis of the Data 

The data are first examined without prejudice as to the origin of any possible excess 
signal by performing a series of statistical tests on each data set. 

The simplest test that can be performed on the data is a calculation of the value of 
X^ for each chop and each field, the results of which are shown in Table 0. All chops have 
values of x^ that lie within the 95% (~ 2a) probability range for 46 degrees of freedom. 
However there appears to be a systematic bias towards x^ values less than 46, indicating 



that the diagonal correlation assumption is incorrect. This is the effect of the residual 
correlation across the co-added data introduced by the presence of atmospheric 1// noise. 
The implications of this correlation are considered further is § |6.3| . 

A maximum likelihood analysis ([Lawrence et al. 1988|) has been carried out on each 
data set to determine the most likely amplitude for any excess variance above the noise 
term described by the uncertainties, cxj. For simplicity, we assume that all points within a 
single data set are independent, ignoring the correlation between adjacent data points that 
is introduced by smearing of any true CMB signal by the beams, or by the presence of 1// 
noise in the data. 

Each data set is analyzed separately using a maximum likelihood estimator with the 
likelihood of an excess variance of al being: 
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The normalized likelihood as a function of (Xe for each difference within each field is shown 
in Figure ^ There is no data set for which o"e 7^ is significantly more likely than o"e = 0, 
as expected from the x^ analysis. To determine limits on ae from these curves, we adopt a 
Bayesian approach with a suitable choice of prior. For consistency with experiments that 
span a similar range of angular scales (e.g. [Myers et al. 1993| ) we adopt a prior that is 
uniform in a^. The 95% ((T95 ~ 2cr) and 99.7% (0-99.7 ~ 3cr) limits are then calculated using 
the highest probability density method (HPD, Perger 1985[ ), in which 



100 X 



CL(ae)dae 



J^ L{ae)dae 



with the constraint that L{ai) = L{au), yields / = 95 or 99.7%. The values ai and au are 
the lower and upper confidence limits respectively. The 95% and 99.7% confidence limits to 
ai are zero for all of the data sets presented here. The 95% and 99.7% confidence limits for 
au are summarized in Table |^. The position of the likelihood peak is also indicated. 

We now examine the validity of the assumption that any excess variance is uncorrelated 
across the data set. The autocorrelation function of the data, Cr, defined as: 



^r — {Vi Vi+r) 

N-r 

E {yilaf){yi+rl(yf+r) 



(9) 
(10) 
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( where yi and cij are defined in equations |^ and ^ with an associated uncertainty, <Jc,r, 
given by: 

4,. = WZT^ (11) 

1=1 
is shown in Figure ^ for data from Field 2. Significant structure can be seen in some of the 
differences. Possible non-astronomical sources for this structure are now considered. 



6. Limits on Systematic Effects in the Data Sets 

6.1. Systematic effects correlated with scan time 

The standard SuZIE observing mode, in which the starting RAs of drift scans are 
alternately 12' and 18' in RA ahead of the nominal source position (see Figure jT]), allows 
checks for systematics that are a function of time within a scan rather than pointing 
location on the sky. To check for such effects, the cross-correlation of the co-added data 
corresponding to the 12' offset with the co-added data corresponding to the 18' offset has 
been calculated for each chop. Artifacts that occur at the same time after the beginning 
of a scan will be seen in the cross-correlation function as a peak at zero time-lag. True 
astronomical signals in the data will yield a feature in the correlation function at a lag of 
At = —AAO cos 6 where A6 = 6' and 6 is the source declination. Since both fields are close 
to the celestial equator At ^ —24 s. 

Figure ^ shows the cross-correlation of data from the two offsets for data sets 
corresponding to Field 2. None of the chops show any strong features at zero lag, indicating 
that there are no systematic effects in the data that occur at the same time after the 
beginning of a scan. Peaks in the correlation functions at At 7^ are seen, but there are 
no peaks at At = —24 s that would indicate the presence of true astronomical signal (the 
peak in t456 is at —30s and is not repeated in any other chop). The most likely cause of 
these peaks is chance correlation of low-frequency 1// noise which, in this experiment, is 
dominated by atmospheric emission fiuctuations (see § |6.3|) . The long (~ 5 s) coherence time 
of atmospheric emission fiuctuations correlates data at the end of one scan with the data 
at the start of the next. The telescope moves by about 30' between scans, corresponding to 
a change in position of the beams of < 8 m at heights < 1 km above the telescope. This 
is less than the correlation length of atmospheric structure which is ~ 25 m (assuming a 
wind speed of 5ms~^ and a measured correlation time of order 5 s). Co-adding many scans 
will reduce the magnitude of this correlation but will not eliminate it entirely. This effect 
manifests itself in the cross-correlation functions shown in Figure ^ as non-zero correlations 
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observed at time lags of ±120 s. 



6.2. Cold stage temperature drifts 

During SuZIE observations, the telescope is stationary during a scan, then re-acquires 
the source before the next scan begins. The motion of the telescope between scans results 
in a small excursion in the temperature of the 300 mK stage at the start of a new scan. 
The output of the temperature sensor on the 300 mK stage, co-added in the same way as 
the science data, is shown in Figure |^. The amplitude of the temperature excursion at the 
beginning of the scan is ~ 80 nK. The temperature recovers with a time constant of ~ 10 s, 
and is stable to ~ 10 nK rms for the duration of the scan. 

Although the excursion at the beginning of each scan is very small, we have carried out 
a number of checks to determine whether this temperature drift causes residual artifacts in 
the co-added data sets. First, as described in § |6.1| , such an artifact would be seen as a peak 
at zero lag in the cross-correlations of data from the two RA offsets, shown in Figure |^. 
No significant peak is seen. Second, since the temperature sensor data is sampled at an 
identical rate to the science data, the measured temperature drift can be correlated with 
the science channels on a scan-by-scan basis. The maximum correlated change in signal in 
the science channels is ATcmb < 4/iK, well below other sources of noise. 

This analysis assumes that the shape and the phase of any temperature change in 
the science channels is identical to that observed in the temperature sensor. Because of 
differences in the thermal properties of the sensor and the bolometers, this may not be 
true. Consequently, there could be a large residual effect in the science data that is poorly 
correlated to the sensor data. As a final check, new co-added data sets were generated after 
excluding the first 21 s (7 bins) of data from each scan, and the likelihood analysis in §^ 
repeated. As shown in Figure P, the sole effect on the calculated likelihood is an increase in 
the confidence limits caused by the effective reduction in integration time. 

From these various tests, we conclude that any signal in the data caused by drifts in 
the temperature of the 300 mK stage is much smaller than the experimental uncertainties 
and can be ignored. 



6.3. Atmospheric effects 

Fluctuations in atmospheric emission, caused predominantly by variations in water 
vapor content, contribute 1// noise to the data with a typical correlation time of 5 s, 
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significantly less than the scan length. We therefore assume that data corresponding to the 
same coordinates on the sky but separated in time by more than one scan are uncorrelated 
(this assumption ignores the correlation of data at the end of one scan with data at the 
beginning of the next that causes correlations on angular scales of 24' and 36' in the full 
co-added data sets; we assume that the effects of atmospheric noise correlated between 
scans are small). This assumption allows an estimate of the residual contribution to the 
correlation function of the co-added data. 

Representing the ith data point in the jth scan by the sum of ti, the true astronomical 
signal that remains unchanged between scans, and riij, the noisy signal from all sources 
for which the coherence time is less than the scan length, the correlation function of the 
co-added signal is then: 

Sr = Cr+K (12) 

Here Cr = {ti tj+r) is the correlation function of the true astronomical signal and A/'r is the 
residual correlation function of the noisy part of the signal. It can be shown (Appendix A) 
that Afr is given by: 

K - ^^^jj-^ (13) 

J 

where A/'rj is the correlation function of riij in the jth scan. This expression can be 
rewritten as: 
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(14) 



where the first term is the correlation function of the atmospheric noise in a single scan 
averaged over all scans, and the second term is equal to 1/A^ if the uncertainty, aj, is 
roughly constant over all scans. Since the correlation function of a single scan is given by: 



SrJ — {[ti ~\~ ni,j)[ti^r + fT^i+rj)) 

\^J- I J\fj 



(15) 
(16) 



it is clearly not possible to calculate Afrj independently of the contribution from the true 
astronomical correlation function. However, if we assume that the atmospheric noise within 
the individual scans dominates the true astronomical signal in each scan then Afrj ~ Srj 
and the estimated residual correlation function from atmospheric noise, N"^, can then be 
calculated using Equation |13| with Afrj = Srj. Note that Mj-j ^ Cr is implicit in the 



assumption made in §||, Equations ^ and |], that the uncertainty, aij, associated with the 
data points in an individual scan is equal to aj, the rms of an entire scan. 
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How valid is this assumption? Since the residual atmospheric noise calculated using 



Equation |1^ is approximately equal to the average correlation function per scan divided by 
the total number of scans, the contribution of the true astronomical correlation function, 
Cr, to M^ calculated using Equation |TB| is ~ Cr/N. Since N is typically 100 per RA offset for 
each field, this contribution is clearly very small compared to Cr- Therefore, if the calculated 
function, J\f^, is compared with the measured correlation function of the co- added data 
there are several possibilities: (i) if Afrj and J\fr ^ Cr then the contribution of Cr to J\f^ will 
be negligible and N'r ^ill be a good fit to the correlation functions shown in Figure §; (ii) 
if Mr J and AC ^ Cr then the calculated function Af^ will be ~ Cr/N. Thus Afr will have 
the right shape, but will have an amplitude that is too low by approximately two orders of 
magnitude; (iii) if A^. and Cr have comparable magnitudes then A/^^ will be a poor fit to the 
measured correlation functions. 

Figure ^ shows the estimated correlation function of the residual atmospheric noise, 
A/r^ calculated using the methods given above, overlaid on the correlation function of the 
co-added data corresponding to one RA offset taken from Field 1. At angles less than 10', 
the correlation function of the co-added data is moderately well fitted by Af^, corresponding 
to case (i) above and indicating that the majority of the observed correlation in the data 
on these scales is due to atmospheric noise or other sources of 1// noise. There are some 
large deviations however, particularly in the data from row 2 (pixels 4, 5 and 6). Since this 
row is known to be more susceptible than row 1 to atmospheric fluctuations, these large 
correlations are likely to be the result of atmospheric drifts that are correlated on time 
scales longer than one scan (see §|6.1|). 



7. Comparison of the Data with Gaussian Models of CMB Anisotropies 

In order to compare the data to specific models for the power spectrum of CMB 
fluctuations, the correlation properties introduced by the beam response on the sky, and 
the overlap on the sky of the single and double differences, must be taken into account. We 
represent the correlation function of true temperature anisotropies by: 

Cr = C(n,nfc) = (T(n,)T(nfc)) (17) 

where n^ and n^ are directions on the sky. If the CMB sky is assumed to be sampled from 
a gaussian random field then C{na n^) can be written as ATq x J-'{6) where the spatial 
dependence of the function J-'{6) depends only the angular separation, 6, between the two 
vectors n^ and n;,. The model correlation function, Cr, of the measured data in the absence 
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of any noise sources is then given by: 

Cr = Cr = {Ti Ti+r) = AT^ X T{rAe) (18) 

where T, and Tj+j. are the predicted thermodynamic temperatures of pixels i and i + r in 
the data stream after convolution of the true temperatures with the beam of the instrument. 
The quantity A^ is the sampling rate on the sky. For instance, the correlation function 
between data points measured in the single difference data set, d^i, is given by: 

{TiTi+r)3i = I /(T(na)T(nfe)) [v43(ni-n„)-v4i(ni-n„)] [A3(ni+,,-n,,)-v4i(ni+,.-nfc)] dn^dri;, 

(19) 
where Aj is the beam response function of pixel j in the SuZIE array. 

A simultaneous likelihood analysis for all data from one field can be performed by 
combining all four differences into one data vector y of length AN, where N is the number 
of 45" data bins. The likelihood of the model correlation function, given the data, can be 
calculated as follows: 

^^^^°' ^^ = (2^)iv/2|M|i/2 ^M-y^M-' y) (20) 

The 4N X AN matrix M is composed of sub-matrices m^^ where p,q = 1, 4. The elements 
of each sub-matrix are given by: 

K]=c^r.+-^-. (21) 



where i,j = 1, N. The function C_j is given by Equation |T9| with the beam responses for 
the appropriate differences inserted. The function A/J^j is the correlation function of the 
noise component to the data. 

We need, however, to account for the loss of degrees of freedom introduced by the 
removal of a mean and linear drift from each scan. The process, adapted from Bond et al. 
(1991), is summarized here. The true CMB signal that has been subtracted from the data 
by the removal of an offset and linear drift can be represented as an unknown amount, a, 
of a set of fitting functions, F. The data vector, y, can be corrected for this process by 
defining a vector ytrue such that: 

y = ytruc - a^F (22) 

The unknown vector a^ has length 8 (the product of the number of differences and two 
fitting functions); F is a known 8 x 4A^ matrix with elements equal to the value of the 
fitting functions {Fij = 1 for i = 1, 3, 5, 7 is the basis function for the offset, and Fij = j 
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inodulo A^ for i = 2,4,6,8 is the basis function for the slope, where j = 1, 4A^). Thus the 
on in Equc 

L{ATo,J^) 



hkehhood function in Equation [2^ should be written as: 

1 



Assuming a uniform prior for, and integrating over, the unknown vector a yields: 



(23) 



Lnew(ATo, J") OC 



(|M||F^M-iF|)V2 



exp 



y^M-iy + y^M-iF(F^M-iF)-iF^(M-i)^y 



(2' 



Extending this analysis to include data from both fields is quite simple since the two 
patches of sky are widely separated on the sky, and also widely separated in terms of 
observing time. Thus the two data sets are completely uncorrelated both spatially and 
temporally and the matrix M for the complete data set is: 



M 



tot 



Ml 
Ms 



The likelihood function for all data from both fields is then: 



^tot 



Li X L2 



(25) 



(26) 



where Li and L2 are calculated from Equation 2^. The calibration uncertainty is included 
at this point in the following manner. The measured amplitude of the sky correlation 
function, ATq is assumed to be related to the true value ATq by ATq = G x ATq, where G 
is gaussian-distributed with width o"g (the fractional calibration uncertainty). It can then 
be shown ( [Ganga et al. 1997b| ) that the likelihood of the true amplitude is obtained from 
the likelihood of the measured amplitude by: 



L(Aro,^) 



1 



2n gq ATq 



X 



d{Nr^) exp 



-(AT^ 



ATn)2 



2a^ATo 



L(Ar;^) (27) 



There remains the issue of the noise term A/J^j in Equation ^ The simplest 
assumption is that the noise terms are uncorrelated between separate pixels and differences 
and that A/J^j = 5{p ~ q,j — i) x (crf)^ where af is the uncertainty associated with the ith 
data point in difference p. However, as shown in § |6.3| , there is a significant component of 1// 
noise in the data that introduces correlations between pixels and also between the different 



chops. In principle, the non-diagonal noise terms can be estimated using Equation [13. 
However, this will not account for the residual features in Figure that are believed to be 



-15- 



caused by slow drifts correlated across more than one scan, and that are not well fitted by 
Equation 0. 

Instead we adopt the following approach. Because N^ is such a good fit to the 
correlation function of the co-added data, there is unlikely to be a high signal/noise 
detection of astronomical signal. Thus the data can be analyzed assuming the simple 
(5-function form for A/J^j with the constraint that the results from the likelihood analysis 
should be treated as upper limits only to any signal in the data. The possible effects of 



ignoring the correlated noise component are discussed in §7.1 



To summarize, the elements of the sub-matrices m^'' are mH- 



j — i and A^ = 0'75 x cos 5 is the sampling interval of the data 
on the sky. The modified correlation function, C {rA6) is obtained from Equation IS. 



^{p ^ Qj f') ("^p)^ where 



7.1. Comparison With a Gaussian Autocorrelation Function 



The generalized form for C{6) is an expansion in terms of Legendre polynomials, 
P^cos^): 

-| oo 

CiO) = —Y.i2l + l)CiPiicose) (28) 



An 



1=2 



where the function Ci depends on the assumed cosmological model. The functional form 
of Ci due to anisotropies imprinted at recombination has been extensively modeled in the 
literature ([White et al. 199^ , [Bunn and Sugiyama 1995| , [Gorski et al. 1995i [Ratra et~aT 



1997) and the implications of the SuZIE measurements for these models are considered in 



a companion paper (|Ganga et al. 1997a]) . In this paper we consider our measurements in 
terms of secondary sources of anisotropy from a population of S-Z clusters, since these may 
be comparable to, or larger than, primary anisotropies at these angular scales ( Pond 1995| ). 
The analytically simple gaussian autocorrelation function (CAGE) model for C{6), while a 
poor approximation to the power spectrum of primordial anisotropies, is a useful tool in 
this case. 



The functional forms for the GAGE is: 



C{9) =Coexp 



2^2 



(29) 



with two free parameters, Cq = ATq, the variance, and Oq, the coherence angle, of the 
fluctuations. In the absence of other noise sources, observations of anisotropies distributed 
with this correlation function, using a gaussian beam of EWHM 9b, yield measurements 
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with correlation function: 



-e'' 



2(2/32 + ^2) 



(30) 



where (3 = 0fe/V81og2. 



The SuZIE beams are assumed to be gaussian, with a FWHM of 1'7, derived by 
normahzation to the measured sohd angle. If a gaussian is fitted directly to the beam 
profiles, then the derived FWHM is 6% larger, however, this gaussian overestimates the 
beam shape at several beam widths away from the beam center. The 6% difference can 
be taken as some measure of the uncertainty that is introduced by assuming a gaussian 
beam. Ganga et al. (1997b) show that beam uncertainties of this magnitude do not affect 
the results of the modeling although they must be included in the calibration uncertainty, 
as we have done in ^. 

The likelihood contours as a function of (Cq)^/^ and coherence angle, ^o^ from a 



simultaneous likelihood analysis of all data from all fields, are shown in Figure 10a. A 
peak in the likelihood function at non-zero values of Co 7^ is obtained for all coherence 
angles less than 10'. Using the HPD method with a uniform prior in (Cq)^^^, the 68%, 
95% and 99.7% confidence limits on {CqY^'^ for each value of 6^0 have been calculated, and 
are shown in Figure |Tn|b. The peak in the likelihood function has low significant (< 68% 
confidence) and is consistent with noise. Consequently we confine ourselves to considering 
the 95% and 99.7% confidence limits. The coherence angle at which maximum sensitivity 
to Co is obtained is 6*0 = I'l; at this coherence angle, (Cq)^''^ < 58 /iK (95% confidence) and 
< 72 /iK (99.7% confidence). Expressing these limits in terms of AT/T yields: 

^^ < / 2.1 X 10-5 (95% confidence) 

' - \ 2.6 X 10-5 (99.7% confidence) ^ ' 

These limits have been calculated without including the correlated terms in the matrix 
M that are introduced by the presence of atmospheric noise in the data. As such we believe 
the limits to be conservative. To check this assumption, we have re-analyzed just the data 
corresponding to the d^i and dQ4 chops, at the coherence angle of maximum sensitivity 
(^0 = I'l) oiily, with the inclusion of the terms N'fli in Equation |2l|. These terms have 



been calculated as outlined in § |6.3| . Ignoring data from the triple beam differences, ^123 and 
^456, does not significantly degrade the derived upper limits. We find that the limits that 
are derived when the correlated noise terms are included are somewhat lower than those 
calculated when the correlated terms are ignored. Consequently the assumption that our 
upper limits are conservative is valid. 
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It is also useful to express these limits in terms of C(0) = Cq6q/{2I3'^ + 6^), the variance 
of sky fluctuations described by a GACF after convolution with a single gaussian beam. 
The most likely value, and the 95% and 99.7% limits, for ATrms = [C'(O)]^''^ are shown in 
Figure |TT|. At small values of 6*0, our experiment is unable to distinguish between GACFs 
with different coherence angles and so ATrms tends towards a constant value as ^o tends to 
zero. At large values of 6*0, the SuZIE differencing scheme decreases the sensitivity to ATrms 
and so the limits rise steeply with increasing ^o- Selecting the value of Oq = I'l (the GACF 
to which this experiment is most sensitive) yields upper limits on ATrms of: 



1.6 X 10"5 (95% confidence) 
1.9 X 10^5 (99.7% confidence) 



Ar,„./T < \ :■: : ::_, )z:::r:r:rL m 



As 9q tends to zero, the limits become: 



r 1.2 X 10-5 (95% confidence) 
™^/ - \ 1.5 X 10-5 (99.7% confidence) ^ ' 



8. Discussion 

These results are compared with other limits on similar angular scales and are used to 
set limits on sources of secondary anisotropy. 



8.1. Comparison with other measurements at arcminute scales 

Several other experiments have explored CMB anisotropies on similar angular scales 
to the SuZIE measurements. In general, the most sensitive limits prior to this work have 
been obtained from measurements at centimeter wavelengths. These include: (i) the 
NCP (North Celestial Pole) and RING measurements made with the Owens Valley Radio 
Observatory (OVRO) 40 m dish at a frequency of 20 GHz ( [Readhead et al. 1989| , [Myers et 



|al. 1993|) ; (ii) a measurement made with the most compact configuration of the Very Large 
Array at 15 GHz ([Fomalont et al. 1993|) ; (iii) a limit set with the Australia Telescope at 



8.7 GHz ( Subramanyan et al. 1993 ). At millimeter wavelengths (90 GHz), an upper limit to 



arcminute scale anisotropy has been obtained by the White Dish experiment ( [lucker et al. 



19931) 



Of these measurements, only the Owens Valley RING measurement has detected 
structure. Myers et al. (1993) report a detection of structure with 2.3 < AT^tas/T < 4.5 x 10-^ 
(95% confidence limits) that is inconsistent with the 95% upper limit of 1.7 x 10-^ obtained 
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with the same instrument in the NCP experiment (the RING measurement contains 96 108" 
FWHM fields compared to the 8 108" fields observed in the NCP experiment). Confusion 
from radio point sources is postulated as the cause of this disagreement. 

The 95% upper limits to a GACF obtained from each experiment, and the SuZIE 95% 
upper limits, are plotted as a function of 6q in Figure |l^. At millimeter wavelengths, the 
SuZIE data sets the most sensitive upper limits to GACFs with 6q < 3'. The SuZIE 95% 
upper limit also lies below the RING detection for coherence angles less than 2'5, consistent 
with the hypothesis that the RING detection is due to a population of radio point sources. 



8.2. Limits on populations of Sunyaev-Zel'dovich sources 

Inverse-Compton scattering of CMB photons by the hot gas in galaxy clusters causes 
the Sunyaev-Zel'dovich effect ( |Sunyaev and Zel'dovich 1972| , 1980). The spectral distortions 



that comprise the S-Z effect can be conveniently sub-divided into a component arising from 
the thermal motions of the electrons in the hot gas, and a component arising from the 
kinematic velocity of the cluster with respect to the CMB rest frame ( [Rephaeli and Lahav 
19911) . 



The non-relativistic form of the thermal effect can be written as a frequency-dependent 

AT/TcMB distortion: 

AT r X 1 

— = - X coth 4 y (34) 

TcMB I 2 J 

where x = hu/kTcuB and the Compton y-parameter is given by y ~ rkTe/rrieC^ for a cluster 
with optical depth r to Compton scattering, and average electron temperature, T^. The 
thermal S-Z effect thus appears as a decrement in the CMB temperature at wavelengths 
longer than ~ 1.4 mm and an increment at shorter wavelengths. Clusters with strong X-ray 
emission typically have values of r = 0.01-0.02 and T^ = 8-15 keV. In the non-relativistic 
limit the cross-over occurs at A = 1.38 mm, otherwise the exact point of zero-crossing 
depends on well-defined relativistic corrections ( [Rephaeli 1995] ). 



The kinematic effect has a spectral signature identical to that of primordial CMB 

fluctuations and is given by: 

AT V 

(35) 



TcMB C 

where v is the peculiar velocity of the cluster with respect to the CMB rest frame. In Cold 
Dark Matter models with density parameter Q = 0.3, the rms value of peculiar velocities 
is predicted to be 400kms~^ with 10% of clusters having i;pec > 700kms~^; for fi = 1, 
the rms velocity is 700kms~^ with a high velocity tail out to 2000 km s~^ ( Pahcall, Cen 
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and Gramann 1994). Therefore, for a bright cluster with r = 0.02 and Tg = lOkeV, the 
kinematic effect is expected to be at least an order of magnitude smaller than the thermal 
effect, except near A = 1.4mm where the thermal effect is zero. 

An important feature of both S-Z components is that the spectral shapes and 
intensities are independent of redshift for clusters with identical properties. Thus theoretical 
predictions of number counts of S-Z sources depend only on models of large-scale structure 
and the evolution of cluster gas properties with redshift. 

Several authors consider the contribution to AT /T anisotropy from a population of 
X-ray clusters. Bond (1995) and Markevitch et al. (1992, 1994) have simulated maps 
of several square degrees of sky, showing the expected contribution from the S-Z effect. 
De Luca, Desert and Puget (1995), Markevitch et al. (1994), Bartlett and Silk (1994) 
and Barbosa et al. (1996) have calculated the expected number counts of S-Z sources 
as a function of flux and have explored the dependence of this statistic on the assumed 
cosmological model. Bond (1995) and De Luca et al. (1995) consider the contribution from 
both the thermal and kinematic effects; Markevitch et al. (1992, 1994), Bartlett and Silk 
(1994) and Barbosa et al. (1996) consider the thermal effect only. 

The correlation function of cluster-induced anisotropics is strongly dependent on 
the assumed cosmological model (see Figure 7 of Pond 19"95| for theoretical models of 
the cluster-induced Ci power spectrum). Experimental results such as those obtained 
with SuZIE are best compared with theoretical models of such anisotropics via detailed 
simulations of model skys, convolved with the appropriate beam response. However, since 
there is no significant detection of structure in the SuZIE data, a simplistic comparison of 
the AT/T limits from the GACF analysis to the models serves as a useful guide to how 
closely current experimental limits are approaching theoretical predictions. 

Since the SuZIE beam size is well matched to the size of a typical cluster at z > 0.1, 
we use the 99.7% confidence limit of AT/T < 2.6 x 10"^ (from Equation ^1]) obtained 
using a GACF with a coherence angle of I'l. Converting to a value for the y-parameter 
yields Ay < 2.5 x 10~^ at the 99.7% confidence level. Ignoring the effects of the SuZIE 
differencing scheme yields a limit to the S-Z flux at 142 GHz from any cluster that may be 
in the SuZIE fields of |A/s_z| < 7.3 mJy. Scaling theoretical S-Z number counts at 400 GHz 
( [Barbosa et al. 1996| ) to 142 GHz yields a prediction of 1.2 clusters with absolute flux 



greater than 7.3 mJy in the 0.06 square degrees observed by SuZIE if fi = 0.3. If fi = 1, 0.2 
clusters per 0.06 square degrees are expected. Thus in a low-fi universe, the SuZIE limits 
are comparable to the theoretical predictions. Because the expected number of clusters is 
small, observations of larger sky patches to a similar sensitivity level are required to allow 
meaningful comparisons with models of the S-Z cluster background. The results reported 
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here indicate that such a test is within reach in the near future. 

9. Conclusions 

Based on observations of two patches of sky covering a total area of 0.06 sq. deg. at 
142 GHz with 1'7 resolution, we have set 95% confidence limits of AT/T < 2.1 x 10"^ for 
CMB anisotropics distributed with a gaussian autocorrelation function with a Ifl coherence 
angle. These limits do not include the effects of atmospheric correlations in the data, but 
a partial re-analysis that accounts for these correlations has been shown to reduce the 
upper limits. Consequently the numbers we present here are conservative. These limits 
are comparable to the best limits obtained from centimeter-wave observations on similar 
angular scales. Because the SZ effect is brighter relative to CMB anisotropy at 142 GHz 
than at centimeter wavelengths, these upper limits on CMB anisotropy provide the most 
sensitive probe of the background of SZ fluctuations to date. 

The upper limits that we have obtained are comparable to fluctuations from the 
Sunyaev-Zel'dovich effect in galaxy clusters predicted for a low-O universe. Observations of 
larger areas of sky at a similar level of sensitivity, coupled to detailed simulations of model 
skies that include the effects of our scan strategy, will significantly constrain models of the 
formation of galaxy clusters. 

An upgraded SuZlE instrument has now been constructed and is currently being 
commissioned at the CSO. This multi-frequency instrument makes simultaneous 
measurements at 142, 217 and 268 GHz and thus allows atmospheric noise, which limits 
the sensitivity of the current system, to be subtracted via correlation between different 
frequency channels. The multi-frequency capability will also allow separation of primary 
CMB anisotropics from secondary fluctuations such as those caused by S-Z clusters. This 
instrument will achieve sensitivity levels significantly better than the results reported 
here, enabling us to detect both primary and secondary anisotropics if they exist at levels 
predicted by current theories. 

This work has been made possible by a grant from the David and Lucile Packard 
Foundation, and by a National Science Foundation grant #AST-95-03226. We thank 
Anthony Schinckel and the entire staff of the CSO for their excellent support during the 
observations. The CSO is operated by the California Institute of Technology under funding 
from the National Science Foundation, Contract 7^AST-93-13929. 
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A. Derivation of the residual atmospheric correlation function 

If the data in the ith pixel of the jth scan is given by yij = tj + riij, where ti is 
the astronomical signal and riij is the noise term, then the co-added signal as given by 
Equations |^ and |] can be written as: 

yi = ti + ni (Al) 

where: 

rii = — (A2) 

3 

We ignore for the moment the simplification that can be applied to this particular data set, 
that (Tjj- = (Tj and is identical for all points within a single scan. 

The correlation function of the co-added data, Sr = {yiVi+r), is then: 

^r = \titi+r) + iP'i^i+r) \^'~') 

= Cr+Mr (A4) 

where C^ is the correlation function of the true astronomical signal and Mr is the residual 
contribution to the correlation function from noise correlated within a scan. Substituting 



the full expression for rii into the definition of the correlation function given in Equation |T0 
yields: 

■^^ - N-r ( \ ( \ ^^^> 

E ( E i/<- j ( E ^l^+rA 

Since we assume that rii^ is uncorrelated between scans, the numerator of this 
expression is zero except when j = j'. Thus: 

N-r , •. , 

Reversing the order of summation in the numerator yields: 

N-r 



EA/;,, X E (l/<)(l/a^.,,) 



K = ^r-7 '"' ^ / r- (A7) 

g f E l/alA f E l/cTlr,A 
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where Mr,j is the correlation function of the noise in the j'th scan and is given by: 

N-r 

A/;. = "^hr^r (A8) 

Since in this data set atj is constant for all values of i, Equation ^^ simplifies to: 

EAr,,,/aJ 
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Table 1. Locations of the two fields observed in April 1994 



Equatorial Galactic 



RA (1950) Dec (1950) 1 



Field 1 10'^21'"49" 4°04'23" 240?01 47?89 
Field 2 le^SO'^l?" 5°56'00" 21?35 33?37 



Table 2. Average magnitudes in /xK s ^ of the linear drift removed from each channel. 





Chop 


Mean drift 


Rms drift 






[mKs-1] 


[^Ks-i] 


Field 1 


tl23 


25 


118 


(190 scans) 


d3l 


16 


45 




i456 


49 


197 




dei 


-6 


41 


Field 2 


tl23 


7 


56 


(277 scans) 


dai 


-39 


103 




^456 


-63 


240 




^64 


-28 


113 
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Table 3. Calculated value of x^ for each difference and each field. The probability of the 

X^ value being exceeded is also shown. 



difference x^ (46 d.o.f.) P{> x^ 



Field 1 


^123 


55.6 


0.157 




^31 


29.1 


0.976 




^456 


31.7 


0.947 




^64 


40.4 


0.705 


Field 2 


^123 


39.3 


0.746 




^31 


35.9 


0.857 




^456 


45.5 


0.491 




^64 


33.5 


0.915 



Table 4. Results of maximum likelihood analysis to determine the magnitude of excess 

variance in each data set. 





Chop 


Peak [^K] 


Upper 


confidence limits, 


On [mK] 


Field 


95% 


99.7% 




Field 1 


^123 


47 


98 


132 






dzi 





79 


120 






^456 


1 


98 


148 






^64 





109 


156 




Field 2 


il23 





52 


77 






^31 





95 


140 






^456 


27 


116 


162 






^64 


1 


94 


140 
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Fig. 1. — An illustration of the SuZIE scan strategy. The upper portion shows the 
arrangement of the six pixels in the focal plane. In the lower portion, the drift scan strategy 
is shown with the location of the field centers given in Table 1 indicated. The opposing 
directions of the line hatching indicate the two RA offsets and the overlapping region is 
cross-hatched. 



-28- 




2000 4000 6000 



10* 




-0.005 



0.005 



(y,,j/f^f)/MK-' 



Fig. 2. — a) The distribution of the scan rms values, 0"^, for all 190 scans obtained on Field 1 
with the (^04 difference, b) The distribution of the quantity yjj/o"^ for a single bin from all 
190 scans obtained on Field 1 with the (i64 difference. The solid line represents the normal 
distribution defined by the rms of Ui^j/a'j showing that this quantity is normally-distributed 
with no extreme values or large wings. 
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Fig. 3. — SuZIE measurements of Field 1. The average uncertainty on each point is ~ 130 /iK 
per bin for ti23, 180 /iK for ^31, 220 /iK for t^^^Q and 200 //K for ^64- The instrument response 
to a point source is also shown in each panel. 
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Fig. 4. — Likelihood values for excess variance, o"e, in the SuZIE data, assuming no correlation 
between data points. The likelihood function for each difference and each field is normalized 
to a peak value of 1. 
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Fig. 5. — The correlation function of the co-added data from Field 2. 
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Fig. 6. — The cross-correlation of co-added scans corresponding to the 12' offset with co- 
added scans corresponding to the 18' offset. These data are taken from Field 2. Systematic 
effects occurring at the same time after the beginning of a scan would yield a peak at At = 0. 
True astronomical signal would be seen as a peak at At = —24 s. These values are indicated 
by the dashed lines. 
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Fig. 7. — Co-added data from the 300 inK stage temperature sensor showing the temperature 
excursion at the start of a scan induced by telescope motions. 
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Fig. 8. — A comparison of the likelihood values for excess variance, a^.^ showing that when 
the first 21 s of data from each scan are excluded from the analysis (dashed line) there is no 
significant change compared to the likelihood function for all the data (continuous line) other 
than that expected from loss of integration time. These likelihood functions correspond to 
the ti23 chop from Field 1. 



-35- 



5x10* 



X 
i 



CO 

o 



-5x10* 



-1x10^ 




10 20 

6/arcmin 



105 



5x10* 



,3. 



-5x10* 



-1x10^ 




10= 



5x10* 



1^ 

a. 



10 20 
6/arcmin 



-5x10* 



-1x10= Ll 
30 




10 20 

6/arcmin 



30 



Fig. 9. — Correlation function of one RA offset of co-added data from Field 1. The continuous 
line is the residual correlation function calculated using Equation f^ . 
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Fig. 10. — a) Likelihood contours, normalized to a peak likelihood of 1, obtained by assuming 
a gaussian autocorrelation function for the distribution of CMB anisotropies. Contour levels 

1 /2 

start at 0.1 (labeled) and increase in steps of 0.1. The most likely value of Cg and ^o is 
marked with a cross, b) The continuous line shows the position of the maximum likelihood 
for each value of ^q, the dashed lines show 63%, 95% and 99.7% confidence limits. 
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Fig. 11. — Limits on AT^ms = [C'(O)]^^^ ^-s a function of coherence angle, 6q. The bold line 
shows the most likely value, the dashed lines show the 95% and 99.7% upper limits. 
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Fig. 12. — SuZIE 95% confidence limits (solid line) to a gaussian autocorrelation function 
for AT/T anisotropies, as a function of ^o- The corresponding 95% confidence limits (long- 
dashed lines) from the Owens Valley Radio Observatory (OVRO) NOP measurement, the 
Very Large Array (VLA), the Australia Telescope Compact Array (ATCA) and the White 
Dish (WD) experiment are also shown (for references, see text). The position of the likelihood 
peak obtained from the OVRO RING experiment is also shown (short-dashed line). 



